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Search for Scalar Bottom Quark Pair Production with the ATLAS Detector 

in pp Collisions at yfs = 7 TeV 



(N 

o 

(N 

a. 



i 

43- 



(N- 

>: 

(N. 

m ■ 
oo- 
cn ; 

(N 1 



X: 



The ATLAS Collaboration 

The results of a search for pair production of the scalar partners of bottom quarks in 2.05 fb _1 of 
pp collisions at ^/s = 7 TeV using the ATLAS experiment are reported. Scalar bottom quarks 
are searched for in events with large missing transverse momentum and two jets in the final state, 
where both jets are identified as originating from a bottom quark. In an 7?-parity conserving minimal 
supersymmetric scenario, assuming that the scalar bottom quark decays exclusively into a bottom 
quark and a neutralino, 95% confidence-level upper limits are obtained in the b x — Xi mass plane 
such that for neutralino masses below 60 GeV scalar bottom masses up to 390 GeV are excluded. 

PACS numbers: 11.30.Pb,12.60.Jv,14.80.Ly 



Supersymmetry (SUSY) [lHSl is a theory that provides 
an extension of the Standard Model (SM) and naturally 
resolves the hierarchy problem by introducing supersym- 
metric partners of the known bosons and fermions. In 
the framework of a generic i?-parity conserving minimal 
supersymmetric extension of the SM, the MSSM, SUSY 
particles are produced in pairs and the lightest super- 
symmetric particle (LSP) is stable, providing a possi- 
ble candidate for dark matter. The scalar partners of 
right-handed and left-handed quarks (squarks) can mix 
to form two mass eigenstates. In a large variety of mod- 
els, the LSP is the lightest neutralino (%?) and the mass of 
the lightest scalar bottom quark eigenstate (sbottom, b±) 
can be significantly lower than the other squark masses. 
Consequently, b 1 could be produced with large cross sec- 
tions at the Large Hadron Collider (LHC). In this Let- 
ter, a search for direct sbottom pair production is pre- 
sented, assuming a SUSY particle mass hierarchy such 
that the sbottom decays exclusively into a &-quark and a 
Xi (t>i — > bxi)- The expected signal is characterized by 
the presence of two jets originating from the hadroniza- 
tion of the 6-quarks and missing transverse momentum 

— its magnitude is referred to as E T llss in the following 

— resulting from the undetected neutralinos. Results of 
searches for direct sbottom production have been pre- 
viously reported by the Tevatron 1(| 11 1 and LEP 12 1 
experiments. The search described here extends these re- 
sults using 2.05 fb _1 of 7 TeV pp collision data recorded 
by the ATLAS experiment at the LHC. 

The ATLAS detector [l3| consists of inner tracking de- 
vices surrounded by a superconducting solenoid, electro- 
magnetic and hadronic calorimeters, and a muon spec- 
trometer with a toroidal magnetic field. The inner detec- 
tor, in combination with the 2 T field from the solenoid, 
provides precision tracking of charged particles for \rj\ < 
2.5 141 ] . It consists of a silicon pixel detector, a silicon 



range \rj\ < 4.9. It is composed of sampling calorimeters 
with either liquid argon (LAr) or scintillating tiles as the 
active media. The muon spectrometer has separate trig- 
ger and high-precision tracking chambers which provide 
muon identification and measurement for < 2.7. 

This search uses data recorded between March and 
August 2011 at the LHC. After the application of 
beam, detector, and data quality requirements, the 
data set corresponds to a total integrated luminosity 
of 2.05±0.08 ft)" 1 [H [13]. The data are selected with 
a three-level trigger system that requires a high trans- 
verse momentum (pr) jet and missing transverse momen- 
tum 17j. Events are required to have a reconstructed 



strip detector, and a straw tube tracker that also provides 
transition radiation measurements for electron identifica- 
tion. The calorimeter system covers the pseudorapidity 



primary vertex associated with five or more tracks, con- 
sistent with the beam spot position. 

Jets are reconstructed offline from three-dimensional 
topological calorimeter energy clusters by using the anti- 
k t jet algorithm [3, ^ with a radius parameter of 0.4. 
The measured jet energy is corrected for inhomogeneities 
and for the non-compensating nature of the calorime- 
ter by using pr™ and ^-dependent correction factors ( (2pj 
and references therein). Only jet candidates with pr > 
20 GeV within \rj\ < 2.8 are retained. During the data- 
taking period, a localized electronics failure in the LAr 
barrel calorimeter created an inactive region in the sec- 
ond and third calorimeter layers (Arj x A<fi ~ 1.4 x 0.2) 
in which on average 30% of the incident jet energy is not 
measured. If either of the two selected leading jets fall 
in this region the event is rejected. The loss in signal ac- 
ceptance is smaller than 10% for the models considered. 

Electron candidates are required to have pt > 20 GeV, 
\r)\ < 2.47, and satisfy the 'medium' selection criteria re- 
ported in [2l|. Muon candidates are required to have 
Pt > 10 GeV, \r]\ < 2.4 and are identified as a match 
between an extrapolated inner detector track and one 
or more track segments in the muon spectrometer. Ad- 
ditional requirements are applied to electron and muon 
candidates when defining lepton (£=e, /i) control regions. 
In this case, electrons must have pt > 25 GeV, pass 



2 



the 'tight' selection criteria in |2l| and be isolated - 
the pt sum of tracks within a cone in the (rj, <f>) plane 
of radius AR — 0.2 around the candidate, £pr, must 
be less than 10% of the electron pt- Muons must have 
Pt > 25 GeV, longitudinal and transverse impact pa- 
rameter within 1 mm and 0.2 mm of the primary vertex, 
respectively, and Spx < 1.8 GeV. Following the object 
reconstruction described above, overlaps between jet can- 
didates and electrons or muons are resolved. Any jet 
within a distance AR = 0.2 of a 'medium' electron can- 
didate is discarded. Any remaining lepton within AR — 
0.4 of a jet is discarded. 

The calculation of £J™ ISS is based on the magnitude 
of the vectorial sum of the pt of the reconstructed jets 
(with pt > 20 GeV and \rj\ < 4.5), leptons — including 
non-isolated muons — and the calorimeter clusters not 
belonging to reconstructed objects [22j |. 

Events must pass basic quality criteria against detector 
noise and non-collision backgrounds and are selected with 
at least one jet with p T > 130 GeV and £^ iss > 130 GeV 
to ensure full trigger efficiency, and one additional jet 
with pt > 50 GeV. Both jets are required to be associ- 
ated to the primary interaction and to originate from a 
6-quark using a tagging algorithm that exploits both im- 
pact parameter and secondary vertex information. Jets 
are tagged for \rj\ < 2.5 and the parameters of the algo- 
rithm are chosen such that a tagging efficiency of 60%, 
10% and <1% is achieved for 6-jets, c-jets, and light fla- 
vor or gluon jets, respectively, in tt events in Monte Carlo 
(MC) simulation 23[. Events with identified electron or 
muon candidates or additional jets with pt > 50 GeV 
are vetoed to exploit the topology of the addressed sig- 
nal. The multijet background contribution with large 
£?p ISS , due to the mismeasurement of the jet energies 
in the calorimeters or to neutrino production in heavy 
quark decays, is suppressed by requiring the smallest az- 
imuthal separation between the E™ lss direction and any 
of the two leading jets, A<fi m i n , to be above 0.4. If a 
third jet with 30 < pt < 50 GeV is found in the event, 
its azimuthal separation with the E^ lss must be above 
0.2, where jets with pt below 30 GeV are not used due 
to the large expected contribution from multiple interac- 
tions. In addition, the ratio of E™ lss to the scalar sum 
of _E™ ISS and the transverse momenta of the two leading 
jets is required to be above 0.25. 

A selection based on the boost-corrected contransverse 



mass [241 ]. wct, is employed to further discriminate sbot- 
tom pair production from SM background processes. For 
two identical decays of heavy particles into two visible 
particles v% and y^, and into invisible particles, the con- 
transverse mass [25[ is defined as ([Et(vi) + Et{v2)} 2 — 
[Pt(»i) — PtC^)] 2 ) 1 / 2 - The boost-corrected contrans- 
verse mass conservatively corrects to account for boosts 
in the transverse plane due to initial state radiation 
and hence protects the expected endpoint in the distri- 
bution. In the case of sbottom pair production with 



b x — > bxi, racT is expected to have an endpoint at 
(m(6 1 ) 2 — m(xi) 2 )/™(&i). To maximize the sensitiv- 
ity across the b x — Xi mass plane, three signal regions 
(SR) are defined as a function of the toct threshold, 
with m CT > 100, 150 and 200 GeV. 

Simulated event samples are used to aid in the de- 
scription of the background and to model the SUSY sig- 
nal. Top quark pair and single top production are sim- 
ulated with MCQNL0 [26J, fixing the top quark mass at 
172.5 GeV, and the next-to-leading-order (NLO) par- 
ton density function (PDF) set CTEQ6.6 Sam- 
ples of W+jets, Z+jets with light and heavy flavor jets, 
and tt with additional 6-jets, tibb, are generated with 
ALPGEN [HI and PDF set CTEQ6L1 The fragmenta- 
tion and hadronization for the ALPGEN and MC@NL0 sam- 
ples are performed with HERWIG j3(j, using JIMMY [31 1 
for the underlying event. Samples of Zti and Wtt are 
generated with MADGRAPH ^ interfaced to PYTHIA [33jj . 
Diboson (WW, WZ, ZZ) samples are generated with 
HERWIG. All background cross sections are normalised to 
the results of higher order calculations 34|. The signal 
samples are generated in the MSSM framework at fixed 
b\ and Xi masses, with BR(6 1 — > &Xi)=100%, using the 
HERWIG++ 30] v2.4.2 Monte Carlo program. The SUSY 
sample yields are normalized to the results of NLO cal- 
culations, as obtained using the PRDSPIN0 [35[ v2.1 pro- 
gram. The CTEQ6 . 6M [11] parameterization of the PDFs 
is used and the renormalization and factorization scales, 
A*r f = are se t to the sbottom mass. NLO cross sec- 
tions vary between 12 pb and 0.05 pb for sbottom masses 
in the range 200-500 GeV. The MC samples are produced 
using parameters tuned as described in 38] and are 
processed through a detector simulation [39| based on 
GEANT4 [i3|. Effects of multiple proton-proton interac- 
tions per bunch crossing are included in the simulation. 

The signal efficiencies vary across the b^xl mass plane. 
As an example, for the SR with toct > 150 GeV the ef- 
ficiencies are found to be between 1% and 6% as the 
sbottom mass increases from 200 GeV to 500 GeV, and 
between 6% and 2% as Am = rn r — m~ decreases. The 

sensitivity for Am below 130 GeV is limited by the selec- 
tion in toct and the trigger-driven requirements on the 
leading jet pt and E T ilss . 

The main SM processes contributing to the back- 
ground are top quark pair and single top production as 
well as associated production of W/Z bosons with heavy 
flavor jets — referred to as Z+hi and VF+hf, respectively. 
In particular, the signal region with toqt > 100 GeV is 
dominated by semi-leptonic ti events as a consequence 
of the toct endpoint at 135 GeV in top pair produc- 
tion, whilst the signal regions with toqt > 150, 200 GeV 
are dominated by the irreducible Z+hf production with 
Z — > vv decay, followed by W+hf production with 
W — > tv. Contributions from multijet, diboson and as- 
sociated production of tt with W, Z or additional &-jets 
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are sub-dominant. Non-collision backgrounds were found 
to be negligible. 

The estimation of the main background processes is 
carried out by defining data control regions where each 
background component is dominant. The background 
estimate in each SR is derived through 'transfer factors' 
equivalent to the ratio of expected event yields in the sig- 
nal and control regions estimated using the MC. The con- 
tributions from top and VF+hf production are estimated 
using a transfer factor from a control region where events 
have exactly one electron or muon, E™ lss > 80 GeV, and 
at least two 6-jets with px > 130 GeV and 50 GeV. The 
transverse mass of the (£, E™ lss ) system is required to be 
between 40 GeV and 100 GeV to select events containing 
W -+ £v. The contribution to this control region from 
other SM processes accounts for less than 10% of the total 
and is estimated from simulation. The Z+hi contribu- 
tion is estimated using a transfer factor from a control 
region where events have two opposite-sign same-flavor 
leptons (£ + £~), at least two 6-jets with pt > 80 GeV 
and 50 GeV, and invariant mass of the two leptons mu 
between 81 GeV and 101 GeV (Z-mass interval). In ad- 
dition, the transverse momentum of the (£ + ,£~~, E™ lss ) 
system is required to be greater than 50 GeV. The contri- 
bution from top quark production in this control region 
accounts for about 50% of the total and is subtracted 
using a side-band estimate in two 40 GeV mass windows 
above and below the Z-mass interval. 

The sub-dominant background contribution from di- 
bosons, Zti, Wtt and tibb is estimated using MC simula- 
tion and increases from 1% to 10% of the total SM predic- 
tion as the selection cut on toct increases from 100 GeV 
to 200 GeV. Finally, the residual multijet background is 
estimated using data. A sample of multijet events with 
large E!f lss is constructed starting from multijet events 
with low i?™ lss and smearing the momenta of jets with 
response functions. This prediction is tuned in a multi- 
jet dominated control region where the requirement on 
A0 m i n is inverted [ 1. 1 j . The contribution is found to be 
less than 5% across the SRs. 

The total systematic uncertainty on the background 
expectations varies from 21% to 44%, increasing with the 
mcT selection applied, and is dominated by the uncer- 
tainty due to finite data statistics in the control regions. 
The next dominant uncertainty on the SM estimates de- 
rives from the residual uncertainties on the theoretical 
modeling of the top background. It varies between 10% 
and 15% depending on the SR and is evaluated using ad- 
ditional MC samples: ACERMC 0, [431] for the impact of 
initial and final state radiation, PYTHIA and POWHEG Q 
for choice of fragmentation model and generator. The 
residual uncertainties on the W+hf and Z+hi theoreti- 
cal modeling account for less than 5% of the total un- 
certainty. Finally, the experimental uncertainties on the 
6-tagging efficiency and jet energy scale and resolution 
are also considered; across the SRs they vary between each point in parameter space is adopted as the nominal 



5% and 8% and between 6% and 9%, respectively. 

For the SUSY signal processes, uncertainties on renor- 
malization and factorization scales and on the PDF affect 
the theoretical cross section. PDF uncertainties are es- 
timated using the CTEQ6 . 6M PDF error eigenvector set 
and are between 7% and 16% depending on the sbottom 
mass. The variation of renormalization and factorization 
scales by a factor of two changes the nominal signal cross 

section, er nom , by ±15%, independently of m- . In the fol- 

°i 

lowing, the cross sections calculated with scale settings 
2x/i and /i/2 are referred to as <7 min and <r max , respec- 
tively. The impact of detector-related uncertainties, such 
as the jet energy scale and ^-tagging, on the signal event 
yields varies between 35% and 45% and is dominated by 
the uncertainties on the 6-tagging efficiency. 

Table U reports the observed number of events and the 
SM predictions before the mcT selection and for each 
SR. Both transfer factor and MC estimates are given. 
The data are in good agreement with the SM background 
expectations within uncertainties in all cases. Figure [1] 
shows the measured mcT and E^ lss distributions before 
mcT selection compared to the SM predictions. MC esti- 
mates are used, rescaled to match the total integral and 
systematic uncertainties predicted by the transfer fac- 
tor estimates for Z+hi and the sum of top and W+hf, 
respectively. For illustrative purposes, the distributions 
expected for the MSSM scenario with sbottom and neu- 
tralino masses of 300 GeV and 100 GeV, respectively, 
are added to the SM predictions. The results are trans- 



mcT 


top, W+bi 


Z+M 


Others 


Total SM 


Data 


GeV 


TF 
(MC) 


TF 
(MC) 


MC+JS 









67 ± 10 
(60 ± 25) 


23 ± 8 
(16 ± 9) 


3.6 ± 1.5 


94 ± 16 
(80 ± 35) 


96 


100 


36 ± 10 
(34 ± 16) 


23 ± 9 
(12 ± 7) 


3.1 ± 1.6 


62 ± 13 
(49 ± 25) 


56 


150 


12 ± 5 
(13 ± 8) 


12 ± 6 
(8.3 ± 4.7) 


2.7 ± 0.9 


27 ± 8 
(24 ± 13) 


28 


200 


3.2 ± 1.6 3.9 ± 3.2 
(4.1 ± 3.4) (2.8 ± 1.5) 


1.0 ± 0.9 


8.1 ± 3.5 
(8.0 ± 4.9) 


10 



TABLE I: Expected and measured number of events for an 
integrated luminosity of 2.05 fb _1 . Z+M and the sum of 
top and Vy+hf are estimated using a transfer factor esti- 
mate (TF). The column labelled as 'Others' reports multijet 
background predictions as estimated with a jet smearing ( JS) 
method, and sub-dominant SM backgrounds estimated from 
MC. For comparison the numbers obtained using MC samples 
are shown in parenthesis. The total systematic uncertainties 
are also reported. 

lated into 95% confidence-level (C.L.) upper limits on 
contributions from new physics using the CL S prescrip- 
tion 45[. The SR with the best expected sensitivity at 
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FIG. 1: Measured m C T (left) and E™ iBB (right) distributions 
before the mcT selection compared to the SM predictions 
(solid line) and SM+MSSM predictions (dashed lines). The 
dashed grey band represents the total systematic uncertain- 
ties. 



of experimental acceptance and efficiency. The observed 
(expected) excluded values of a eS at 95% C.L. are 13.4 fb, 
9.6 fb and 5.6 fb (15.2 fb, 9.2 fb and 4.7 fb), respectively 
for m CT >100, 150, 200 GcV. 

In summary, we report results of a search for sbottom 
pair production in pp collisions at ^fs — 7 TeV, based 
on 2.05 fb -1 of ATLAS data. The events are selected 
with large _E™ 1SS and two jets consistent with originating 
from 6-quarks in the final state. The results are in agree- 
ment with SM predictions for backgrounds and translate 
into 95% C.L. upper limits on sbottom and neutralino 
masses in a given MSSM scenario for which the exclusive 
decay b 1 — > bxl is assumed. For neutralino masses below 
60 GeV, sbottom masses up to 390 GeV are excluded, 
significantly extending previous results. 
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FIG. 2: Expected and observed exclusion limits, as well as 
±1(7 variation on the expected limit, in the b 1 — Xi mass 
plane. The band around the observed limit delimited by the 
two dashed lines shows the effect of renormalization and fac- 
torization scale variation. The reference point indicated on 
the plane corresponds to the MSSM scenario with sbottom 
and neutralino masses of 300 GeV and 100 GeV, respectively. 
Results are compared to previous exclusion limits from Teva- 
tron experiments. Results from LEP cover the region with 
sbottom mass below 100 GeV. 



result. Figure [2] shows the observed and expected exclu- 
sion limits at 95% C.L. in the b l — Xi mass plane, as- 
suming BR(6 1 — > 6xi)=100%. Systematic uncertainties 
are treated as nuisance parameters and their correlations 
are taken into account. For the MSSM scenarios consid- 
ered, the upper limit at 95% C.L. on the sbottom masses 
obtained in the most conservative hypothesis, cr min , is 
390 GeV for m~ = 0. The limit becomes 405 GeV for 
cr nom and 420 GeV for er max . Neutralino masses of 120 
GeV are excluded for 275 < mg < 350 GeV. The three 
signal regions are used to set limits on the effective cross 
section of new physics models, er off , including the effects 
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